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ABSTRACT 


The  paper  gives  a  brief  overview  of  work  on  atmospheric  turbulence  at  DFVLR,  and  describes  the  research  tools, 
which  include  instrumented  aircraft  and  numerical  models.  Some  results  of  research  on  turbulence  characteristics  are  given; 
these  are  mostly  in  the  convective  boundary  layer.  The  report  discusses  their  application  to  the  study  of  aircraft  response  and 
airframe  loadings. 

Paper  presented  to  the  Sub-Committee  on  Right  of  Flexible  Aircraft  in  Turbulence,  under  the  chairmanship  of 
Dr  G.Coupry. 


Ce  rapport  fournit  une  synthese  des  travaux  effeciues  au  DFVLR  sur  la  turbulence  atmospherique,  avec  unc 
description  des  moyens  de  recherche  employes,  qui  comprennent  des  aeronefs  equipes  d'instruments  et  des  modeles 
numeriques.  Quelques  resultats  de  recherche  sur  les  caracteristiques  de  la  turbulence  y  sont  indiques;  ceux-ci  concernent 
principaiement  la  couche  limite  convective.  Le  rapport  examine  les  possibility  des  applications  de  ces  resultats  a  1‘etude  de 
la  reponse  des  aeronefs  et  de  leurs  structures  aux  contraintes  de  la  turbulence.  Le  rapport  a  ete  presente  au  sous-comite  pour 
le  vol  des  aeronefs  a  structure  non-rigjde  en  milieu  turbulent,  preside  par  le  Dr  G.Coupry. 
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A  brief  overview  of  ataoepheric  turbulence  work  at  CfVIR  is  presented.  The  main  research  tools 
covering  aircraft  and  numerical  models  are  described.  Sene  results  on  turbulence  characteristics  mostly 
in  the  convective  boundary  layer  are  given,  and  their  use  for  aircraft  response  and  loads  studies  is 
illustrated. 


i.  DfncoDcncM 

This  report  gives  a  brief  overview  cn  the  status  of  atmospheric  turbulence  research  mainly  at  the 
Institute  of  Atmospheric  Physics  of  the  German  Aerospace  Research  establishment  (CffVIR)  -  with  sene 
limited  reference  to  work  going  on  elsewhere.  The  knowledge  of  atmospheric  turbulence,  in  particular  in 
lower  altitudes,  has  advanced  considerably  due  to  berth  observational  and  modeling  studies.  Fran  the 
available  experimental  and  numerical  data  a  variety  of  simple  models  (analytical  formulae)  have  been 
derived.  A  recent  and  comprehensive  review  of  the  state  of  the  art  with  particular  reference  to 
engineering  applications  is  given  by  Fanofsky  and  Dutton  (1984) . 

The  material  presented  may  contribute  to  aircraft  response  and  loads  studies  in  two  ways:  firstly, 
to  (apply  additional  data  sets  from  various  altitude  ranges  and  meteorological  coi'ditians,  and 
secondly,  to  provide  a  basis  consisting  of  a  theoretical  framework  and  a  complete  data  base  far  casting 
different  approaches  (and  possibly,  to  derive  a  unified  procedure) . 

The  paper  consists  of  two  main  parts:  one  to  introduce  briefly  the  various  research  tools 
available  for  atmospheric  turbulence  studies  at  DFVLR,  and  another  to  describe  several  results  that  are 
believed  to  be  of  importance  to  aircraft  response  and  loads  studies. 


2.  RESEARCH  TOOLS 


The  complex  nature  of  atmospheric  turbulence  (see  Figure  1)  makes  it  necessary  to  use  a  variety  of 
research  tools  in  a  complementary  and  synergistic  approach.  The  Institute  of  Atmospheric  Biysics  uses 
several  research  aircraft  and  sene  numerical  models.  Figure  2  shows  the  different  kinds  of  results  that 
can  be  obtained  from  these  tools. 


2.1  Atmoepher  2  research  aircraft 

The  Institute  of  Atmospheric  Ihysics  is  currently  using  eight  research  aircraft  for  various 
meteorological  programs.  Among  these  a  two  engine  jet  aircraft  Dassault  Falcon  E  and  three  instrumented 
ASK  16  powered  gliders  are  fully  (Falcon)  or  partly  (motor  gliders)  equipped  for  the  measurement  of 
atmospheric  turbulence  and  wind. 

Mean 

The  two-engine  jet  is  a  hi^i  altitude  fast  research  aircraft  instrumented  for  wind,  turbulence  and 
cloud  physics  measurements  and  able  to  carry  remote  sensing  equipment.  A  description  of  the  aircraft 
and  some  results  of  measurements  are  given  by  Finpel  (1987) .  A  detailed  manual  was  capiled  by 
Meischner  (1985) .  Figure  3  shows  a  side  view  with  the  locations  of  the  sensors.  Table  1  gives  setae 
aircraft  parameters.  The  following  modifications  have  been  made  to  the  basic  aircraft: 

-  a  nosebocm  of  1.8  m  length  for  the  installation  of  sensors  in  the  least  disturbed  flow  area, 

-  an  attachment  point  for  external  loads  cn  the  bottom  of  the  fuselage,  which  can  be  used  for  travel 
pods  but  is  now  mainly  used  for  FMS  droplet  spectroraeUsrs, 

-  two  photographic  windows  equipped  with  optical  glass,  520  mm  in  diameter  (alternatively,  pressure 
proof  pods  allow  operation  of  sensors  without  windows) , 

-  four  pests,  80  nn  in  diameter,  in  the  top  of  the  fuselage,  thich  are  used  as  inlets 
for  air  chemistry  instruments  and  pyrradiemeters, 

-  a  sealable  aperture  250  am  x  570  an  on  the  left  hand  side  of  the  fuselage  (special  windows  for 
optical  and  microwave  remote  sensors) . 

Part  of  the  instrumentation  far  meteorological  measurements  is  installed  permanently,  another  part 
is  optional  depending  on  the  purpose  of  the  specific  mission.  Table  2  gives  an  overview  of  the 
instrumentation.  Resolution  and  absolute  accuracy  of  th*  parameters  are  given  in  rable 

The  data  are  recorded  by  a  fixed-wired  PCM  unit.  160  channels  are  sampled  with  10  Hz,  16  channels 
(for  turbulence  aBasuments)  are  sampled  with  100  Hz  each  with  a  resolution  of  12  bit.  Measurement 
speed  far  turbulence  studies  is  around  100  tb/b.  This  gives  a  spatial  resolution  along  the  flight  track 
of  1  a  .'or  turbulence  variables  and  of  10  m  for  others. 
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Motor  gliders 

Dub  to  their  limited  altitude  range  and  payload  the  motor  gliders  are  primarily  used  in  the 
atmospheric  boundary  layer.  A  low  flicjrt  speed  enabling  the  resolution  of  very  small  scales  on  cne  hand 
arri  a  high  degree  of  flexibility  on  the  other  hand  makes  then  the  ideal  research  platform  even  for 
complex  terrain  stxxiies.  The  investigation  of  the  atmospheric  boundary  layer  requires  high  spacetime 
resolution,  and  this  is  achieved  by  the  utilization  of  three  identically  equipped  aircraft. 

The  ASK  16  is  a  low-wing  aircraft  manufactured  of  wood  and  metal.  Figure  4  illustrates  the  form 
and  dimensions  of  the  aircraft.  The  aircraft  has  two  seats,  the  second  seat  may  also  serve  to  install 
additional  equipment.  Sane  aircraft  parameters  are  given  in  Table  1. 

All  three  aircraft  pmnrnn  identical  instrumentation.  A  complete  description  of  the  measurement 
system  and  of  data  available  firm  several  field  experiments  is  given  by  Jochum  et  al.  (1984),  a  short 
survey  by  Jochtsn  et  al.  (1987).  Table  4  provides  a  list  of  equipment.  The  sensors  for  temperature, 
humidity  and  pressure  are  imtalled  in  arel  on  the  wing  instrument  pod.  A  pitot  static  tube  protrudes 
forward  from  the  point  of  the  pod  so  that  the  total  (static)  pressure  aperture(s)  are  1440  (1311)  nm 
forward  of  the  wing  leading  edge.  The  two  accelerometers  are  installed  in  close  vicinity  of  the  center 
of  gravity.  The  two  gyros  are  mounted  in  the  cockpit  panel  and  are  equipped  with  pick-offs  for  data 
registration. 

The  speed  range  is  between  75  kny'h  and  170  knyh.  Measurement  flights  are  conducted  at 
approximately  120  km/h.  There,  the  linear  approximation  of  the  aircraft  lift  equation  used  in  the 
vertical  wind  calculation  is  valid. 

The  determination  of  vertical  windspeed  is  based  on  the  aerodynamic  method  described  by  Lenschow 
(1976) .  The  parameters  needed  in  the  aerodynamic  lift  equation  for  the  powered  gliders  ASK  16  are  given 
by  Hacker  (1982).  The  method  was  extensively  tested  by  flying  special  patterns  (pilot-induced  pitch  and 
roll  oscillations)  in  calm  air  with  negligible  vertical  motion.  The  results  (Hacker,  1982)  show  that 
the  manoeuvre-induced  perturbations  are  sufficiently  damped  by  the  method.  Table  5  shows  an  overview 
over  accuracies,  resolution  and  scales  attained  under  standard  measurement  flight  conditions 
(120  kn/h  speed,  20-30  km  flight  legs). 

A  small  digitalization  system  (MINIDIG)  was  developed  at  the  Institute  of  Atmospheric  Riysics. 
There  are  16  channels  available  with  a  resolution  of  12  bits  each.  The  sanpling  rate  can  be  selected  at 
.1,  1,  10  or  100  Hz.  In  general,  data  are  recorded  at  10  Hz  (which  is  consistent  with  sensor  response 
times) .  Thus,  a  spatial  resolution  along  the  flight  path  of  3-4  m  is  obtaiiwd. 

Data 


In  general,  the  time  series  of  basic  physical  variables  are  low-pass  filtered  with  a  cut-off 
frequency  of  2  Hz  (motor  gliders)  or  0.6  Hz  (Falcon)  (about  17  m)  in  order  to  eliminate  high- frequency 
noise. 

The  data  from  the  different  aircraft  are  combined  to  yield  vertical  profiles  of  mean  values, 
variances  and  spectral  characteristics. 

Various  methods  have  been  employed  to  assess  and  maintain  a  hic^i  overall  quality  of  data.  For 
details  see  the  reports  cited  above.  Willeke  (1985)  has  analyzed  interocnpariscn  flights  of  ASK  16  and 
Falcon  and  finds  good  agreement. 

2.2  Nwrioal  models 

Numerical  models  of  different  ocnplexity  have  been  developed  at  the  Institute  of  Atmospheric 
Ihysics  in  order  to  simulate  turbulence  in  the  atmospheric  boundary  layer.  A  one-dimensional  (ID)  model 
with  third-order  closure  (Finger  and  Schmidt,  1986)  is  able  to  predict  the  statistical  behavior  of  a 
homogeneous  convective  boundary  layer.  It  has  prognostic  equations  for  mean  quantities  and  for 
variances  and  covariances  (second  moments) .  A  three-dimensional  (3D)  model  computes  directly  the  motion 
of  the  large  eddies  and  uses  various  parameterizaticns  of  the  sub-grid  scale  turbulence  (Schumann  et 
al. ,  1987) .  Both  models  have  been  successfully  ooepared  to  laboratory  measurements  (in  a  watertank)  arri 
to  aircraft  observations  (Jochum  and  Schmidt,  1987) . 

The  development  of  a  3D  direct  simulation  model  (resolving  even  the  small  turbulence  scales)  for 
stable  homogeneous  shear  flews  (Gera,  1987)  serves  mainly  basic  research  purposes. 


There  are  various  ways  to  describe  turbulence.  The  two  principally  different  approaches  are 
illustrated  in  Figure  2: 

one  realization  of  the  complete  three-dimensional  flow  field 
-  statistical  methods  (involving  seme  kind  of  averaging)  are  used  to  'fctain  information  about  the 

most  probable  behavior: 

-  spectral  analysis  allows  determination  of  characteristic  energy  containing  wavelengths,  of 
turbulence  intensity,  and  of  dissipation  rates, 

-  probability  distributions  give  additional  information  about  skewness  (third  moment  inlicatiig 
the  relation  between  vigorous  updrafts  and  weak  downdrafts)  and  kurtosis  (forth  moment 
characterizing  the  difference  between  peak  gusts  and  mean  variance) . 


Bcaiples  of  results  obtained  from  these  different  approaches  are  given  consecutively.  They  all 
refer  to  the  convective  boundary  layer.  There  the  dominant  contribution  to  turbulence  energy  nreww  from 
the  vertical  velocity  ccnponent,  which  therefore  has  been  chosen  for  presentation.  The  heic^t  of  the 
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ccnvective  boundary  layer  is  defined  as  the  level  of  the  lowest  temperature  inversion  and  is  denoted  by 

h ■ 

HigBft-diiBBBicnal  flow 

A  three-dimensiaral  numerical  model  (Schumann  et  2d.,  1987)  has  been  used  to  simlate  the  details 
of  turbulent  convection.  The  ramerical  grid  size  is  reduced  to  the  scale  of  turbulent  eddies  in  the 
inertial  subrange.  The  bigger  portion  of  the  turbulent  kinetic  energy  is  resolved,  whereas  the  subgrid 
scale  turbulence  is  modelled  by  second  order  closure.  This  technique  is  termed  large  eddy  simulation 
(IBS) .  In  order  to  validate  the  method,  model  results  have  been  compared  to  data  obtained  from  a 
laboratory  experiment.  Figure  5-a  shows  a  vertical  x-z  cross  section  of  vertical  velocity.  The 
convectivw  boundary  lay**r  is  characterized  by  the  domain  with  some  strong  updrafts  and  weaker 
downdrafts,  whereas  the  velocities  are  generally  snail  in  the  stable  layer  aloft.  Figure  5-b  shows  a 
corresponding  horizontal  x-y  cross  section  through  the  middle  of  the  boundary  layer,  which  confirms  the 
picture  of  a  few  strong  updrafts  and  many  weak  downdrafts. 


Spectra 

From  the  LES  results  vertical  velocity  spectra  have  been  computed  at  various  heights.  Figure  6 
shews  an  exanple  for  z/z,  =0.5.  There  the  inertial  subrange  exhibits  the  characteristic  -5/3  slope 
(von  Kannan  spectrum) .  The  conputed  energy  containing  wavelength  at  about  Xm  =  1.3  z.  is  in  good 
agreement  with  the  observations. 

Seme  authors  have  derived  sirple  models  (mostly  analytic  expressions)  for  spectra  obtained  from 
extensive  field  experiments.  Kaimal  et  al.  (1976)  for  exanple  give  universal  curves  for  velocity 
spectra  in  terms  of  normalized  coordinates.  Figure  7  shows  that  vertical  velocity  spectra  for  different 
normalized  height  ranges  z/z,  collapse  into  a  single  curve  at  the  high  frequency  end,  whereas  the 
normalized  peak  wavelength  x  (low  frequencies)  depends  on  normalized  height.  This  height  dependence 
is  further  illustrated  in  Figire  8,  where  the  dashed  line  represents  the  relationship  (Caughey  and 
Palmer,  1979) 

-  1.8(1  -  -  .0003  exp(-^-))).  (1> 


A  further  analytic  expression  for  vertical  velocity  spectra  in  the  lower  half  of  the  ccnvective  and 
neutral  boundary  layer  is  given  by  Ho j strip  (1982) .  Figure  9  shows  the  fanmila  and  the  agreement  with 
measured  data. 

Probability  distributions 

An  exanple  of  measured  vertical  velocity  probability  distributions  at  different  height  levels  is 
given  in  Figure  10  (Jochun,  1985) .  There  it  is  evident  that  the  principal  characteristics  (variance, 
skewness,  kurtosis)  are  height  dependent  and  in  general  ncn-Gaussian.  The  complete  flow  fields  shown 
partly  in  Figure  5-a  and  Figure  5-b  also  clearly  indicate  the  noo-Gaussian  nature  of  the  vertical 
velocity  distribution  far  convective  turbulence.  The  third  moment  w'J  (which  is  a  measure  of  skewnessl 
as  shewn  in  Figure  11  as  a  function  of  normalized  height  attains  a  maximum  value  of  about  .25  w#3 
around  the  middle  of  the  boundary  layer  (where  the  peak  vravelength  has  a  maximum,  too) . 

vertical  velocity  variance  as  a  function  of  height  is  shown  in  Figure  12  (conputed  from  two 
different  numerical  models)  and  Figure  13  (from  observations) .  The  analytical  expression  derived  from 
the  observations  (Wilczak  and  Ihillips,  1986) 

^  =  {z.s(f  rn  -  .aif  ]}•*  (2) 

is  given  by  the  solid  line  in  Figure  13.  The  ccnvective  scaling  velocity  w#  can  be  estisated  from 
routine  meteorological  data  (van  Ulden  and  Holtslag,  1985) . 

Figure  19  shows  the  good  agreement  between  vertical  velocity  variance  profiles  obtained  from  3D  IBS 
simulations  and  motorglider  measurements . 

4.  APPLICATIONS  TO  AIRCRAFT  RESPONSE 

During  all  measurement  flights  vertical  and  horizontal  accelerations  of  the  aircraft  have  been 
recorded  as  well.  Even  without  knowing  the  aircraft's  response  characteristics  (from  a  system 
identification)  or  the  pilot  induced  movements  (which  are  not  yet  measured) ,  acme  conclusions  about  the 
dynamic  reaction  of  the  aircraft  in  turbulent  air  can  be  drawn.  Spectral  analysis  shows  clearly  that 
the  dominant  turbulence  wavelength  has  virtually  no  influence  on  the  dominant  response  frequency  which 
is  around  0.3-0. 6  Hz  at  all  heights  (Figure  16)  and  with  dominance  of  shear  or  buoyancy  production 
(Reinhardt,  1985) .  Inspection  of  probability  distributions  for  the  same  datasets  reveals  marked 
differences  with  changing  height  (Figure  17)  or  turbulence  production  mechanism  (Figure  15)  far 
vertical  wind  velocity  but  virtually  no  differences  for  vertical  aircraft  acceleration. 

The  results  of  a  different  method  of  analysis  are  shown  in  Figure  18.  There  the  contributions  to 
the  total  variance  resulting  from  different  scale  (or  frequency)  ranges  are  conputed  by  appropriately 
filtering  the  time  series  (Jochum  et  al.,  1987).  small  scales  represent  turbulence  of  stochastic 
character,  intermediate  scale  sizes  are  organized  structures  up  to  the  boundary  layer  depth  scale,  and 
large  scales  refer  tp  eddy  sizes  larger  than  the  boundary  layer  depth.  Vertical  velocity  variance  is 
produced  mainly  by  intermediate  size  scales,  whereas  vertical  aircraft  acceleration  variance  is 
produced  essentially  by  small  scales. 
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5.  GOMCUBiaC  AMD  OUEUKK 

A  brief  overview  of  atmospheric  turbulence  work  at  DFVLR  has  been  presented.  This  work 
concentrates  on  low  altitixie  turbulence  in  the  atmospheric  boundary  layer.  There,  the  turbulence 
characteristics  over  homogeneous  terrain  are  fairly  well  understood,  and  can  be  described  in  terms  of  a 
few  parameters  which  can  be  estimated  from  routine  meteorological  data.  The  accuracy  of  this 
descripiticn  is  satisfactory  for  many  application  purposes. 

Further  research  at  EffVLR  is  aiming  at  understanding  the  atmospheric  boundary  layer  over  conplex 
terrain  and  in  chiming  meteorological  conditions.  Model  development  is  continuing  and  will  finally 
lead  to  a  3D  US  model  with  orography  and  clouds.  The  three  motor  gliders  are  presently  ccnpleted  by  an 
air  motion  sensing  and  high-accuracy  satellite  positioning  system.  The  Falcon  jet  will  have  a  new 
avionics  and  inertial  reference  system  soon  which  is  supposed  to  inprcve  wind  art  turbulence 
measurements. 


Both  aircraft  are  presently  being  equipped  with  sensors  measuring  the  pilot-induced  novenents,  and 
a  system  identification  is  being  performed.  Then  the  ccnplete  dynamic  response  of  the  research  aircraft 
in  turixilent  air  can  be  measured.  Analoguosly,  models  of  the  dynamic  response  are  going  to  be  added  to 
the  existing  numerical  models,  so  that  the  aircraft  response  can  be  confuted,  and  results  of 
computations  can  be  compared  to  measured  data. 
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ASK  10 

Falcon 

maximum  take-off  weight  (MTOW) 

765  kp 

13  000  kp 

permanently-installed  equipment 

70  kp 

included  below 

allowance  for  crew  and  fuel 

200  kp 

5  090  kp 

range 

700  Km 

1220  -  1680  km 

(Jepends  on  altitude) 

service  ceiling 

8500  ft  msl 

42  000  ft 

maximum  rate  of  climb 

2  m/s  at  100  km/h  10  m/s 

Table  1.  Selected  aircraft  specifications 


6 


permanently  installed  sensors  are: 

•  Pressure  and  wind: 

-  Flow  angle  sensor  Rosemount  858  J  with 
transducer  Rosemrunt  1201  for  absolute 
(static)  pressure  and  Rosemount  1221  for 
differential  (impact,  angle  of  attack  and 
sideslip)  pressures. 

-  Inertial  navigation  sytem  Litton  72. 

•  Temperature: 

-  Two  probes  Rosemount  102  (fast  and 
slow). 

•  Humiditiy; 

-  Dew  point  sensor  General  Eastern  101 1 . 

-  Relative  Humidity  sensor  Vais^la  Humlcap. 
—  Absolute  Humiditiy  sensor  ERC  Lyman  a. 
The  Humicap  and  the  Lyman  a  sensors  are 
mounted  inside  the  fuselage  in  a  tube.  On  its 
inlet  the  housing  the  Rosemount  temperature 
sensor  is  fixed.  Close  to  the  humidity  sensors 
the  temperature  and  the  pressure  are  meas¬ 
ured  in  the  tube.  So  the  influences  of  the  heat¬ 
ing  of  the  air  and  the  rise  of  the  pressure  in  the 
tube  of  the  humidity  measurements  can  be  cor¬ 
rected. 

•  Liquid  water  content: 

—  Cloud  technology  (Johnson-Williams)  Sen¬ 
sor  CT-10. 

Optional  sensors  are: 

•  Radiation: 

-  Pyranometer  Eppley  PSP  (upward  and 
downward). 

-  Pyrgeometer  Eppley  PIR  (upward  and 
downward) 

-  Radiometer  Barnes  PRT6  (downward). 

•  Cloud  Physics: 

-  Droplet  spectrometer  PMS  FSSP-100. 

-  Droplet  spectrometer  PMS  OAP-230X. 

-  Liquid  water  content  sensor  PMS  CSIRO- 
King. 

•  Aaerosol  Lidar: 

-  Airborne  Lidar  System  DFVLR  ALEX-F1. 


Table  2.  Instrumentation  Falcon  (from  Fimpel,  1987) 


Temperatures 

•  reverse  flow  temperature  probe  (Pt  100  manu¬ 
factured  in-house) 

•  Pt  100  in  humidity  channel  (manufactured  in- 
house) 

•  Barnes  PRT-5  radiation  thermometer  (one  air¬ 
craft  only) 

Humidity 

•  ERC  Lyman-alpha  hygrometer 

•  Vaisdla  humicap  relative  humidity  sensor 

Pressure 

•  Rosemount  static  and  dynamic  transducers  for 
pitot  static  tube 

Aircraft  motion 

•  Sundstrand  accelerometer  (2) 

•  SFENA  vertical  gyro  (pitch  and  roll  angle) 

•  AIM  directional  gyro  (heading) 

Video  camera(s)  looking  forward  and/or  downward 

(optional) 


Table  4.  Instrumentation  ASK  16 

(from  Jochum  et  al.,1987) 


relative  accuracy  of  derived  quantities 

potential  temperature 

±  07  K 

specific  humidity 

±  05  g/kg 

vertical  wind  speed  *) 

+  2  m/s 

height 

±  *3  m  (abs  ) 

*)  includes  only  errors  of  input  parameters 

horizontal  resolution 

3  3  m 

scales 

16  5  m  -  3  3  km 

altitude  range 

30  m  agl  -  2.5  km  msl 

Table  5. 


absolute  accuracy  of  main  parameters 

Parameter 

Resolution 

Accuracy 

Static  pressure 

±0.12 

±1.5  hPa 

impact  pressure 

+0.06 

±1.2  hPa 

Temperature 

+0.02 

±0.5  K 

Relative  humidity 

±0.02 

±2.0  % 

Absolute  humidity 

±0.01 

±0.7  g/m3 

Horizontal  wind  components  *) 

±5.0  m/s 

Vertical  wind  component 

•) 

±0.8  m/s 

*)  Computed  from  several  parameters 

horizontal  resolution 

1  m 

scales 

5.5  m  -  5  km 

altitude  range 

30  m  agl  -  12  km  msl 

Standard  data  characteristics  ASK  16 
(from  Jochum  el  al..  1987) 


Table  3.  Standard  data  characteristics  Falcon 
(from  Meischner,  1985) 


Figure  1.  Turbulence  structure  ot  the  atmospheric  boundary  layer:  schematic  view 
(Irom  Wyngaard,  1986)  ol  (a)  the  convective  boundary  layer  capped  by  an 
inversion  and  (b)  the  stably  stratified  nocturna'  boundary  layer 


r  e  s  u  Its  < - ^  tools 


Figure  2.  Tools  In  atmospheric  turbulence  research. 


tn  rscw  is 


1 :  Tempo ratur 
2:  Humidity 

3:  Johnson-Wllllams  LWC 
4:  Outlet  of  the 
'humidity  tube' 

5:  Pyrgeometer  downward 
6:  Pyranometer  downward 
7:  Static  pressure 


8:  Imped  pressure 
9:  Flow  angle  sensor 
10:  Photo  windows 
(Lldar  system) 

1 1 :  Pyranometer  upward 
12:  Pyrgeometer  upward 
13:  Droplet  spectrometer 
and  CSIRO-KIng  LWC 


Figure  3.  Falcon  with  Instrumentation  (from  Fimpel,  1987). 
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Figure  4.  ASK  IS  with  Instrumentation  (from  Jochum  et  al. ,  1987). 


Figure  5.  Vertical  velocity  component  w.  (a)  in  an  x-z  cross  section  (from  Schumann 
et  al.,  1987);  isoline  increment  0.5w;  solid  (dashed)  lines  represent  positive 
(negative)  values,  (b)  x-y  cross  section  at  0.62 !z,;  solid  (dashed)  lines  denote 
updrafts  (downdrafts);  isoline  increment  .0Q2m/s.  (Schmidt,  1987,  personal 
communication). 
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Figure  6.  Vertical  velocity  spectrum  from  LES  simulation  at  0.6z,  (from  Schumann  et 
al.,  1987);  solid  (dashed)  lines  denote  computed  (measured  In  a  watertank) 
values. 


f,  »  nh/U 


Figure  7.  Normalized  velocity  spectra  from  observations 

expressed  in  mixed-layer  similarity  coordinates 
(Kaimal  et  al..  1976). 


Figure  8.  Normalized  peak  wavelength  from  observations 

(from  Caughey  and  Palmer.  1979).  The  dashed 
line  represents  the  relationship  (1),  h  denotes 
the  inversion  height. 


Figure  9.  Vertical  velocity  spectra  and  variance:  Hojstrup  formula  (from  Hojstrup. 

1982).  The  surface  layer  friction  velocity  u  and  the  Monin-Obuchow  length 
L  can  be  estimated  from  routine  meteorological  parameters  (van  Ulden  and 
Holtsiag.  1985) 
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Figure  10.  Vertical  velocity  probability  distributions  from  observations  at  different 
heights.  The  abscissa  is  labelled  vertical  wind.  m/s. 


Figure  11.  Vertical  velocity  skewness  as  a  function  of  Figure  12. 

height  (from  Schumann  et  at.,  1987).  solid 
(dashed)  lines  denote  computed  (measured  in 
a  watertank)  values. 


Vertical  velocity  variance  from  model  simu¬ 
lations  (from  Finger  and  Schmidt,  1986),  nor¬ 
malized  with  w .  DW  85  denotes  measurements 
in  a  watertank.  Ml  and  M2  are  models  with  dif¬ 
ferent  turbulence  parameterization. 


Figure  13.  Vertical  velocity  variance  from  observations 

(from  Wilczak  and  Phillips.  1986).  The  solid  line 
is  given  by  relationship  (2),  the  heavy  dashed 
line  represents  observations. 
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Figure  14.  Measured  spectra  of  vertical  wind  (left)  and 

vertical  acceleration  of  aircraft  (right):  shear 
produced  (below)  versus  buoyancy  produced 
(above)  turbulence. 


Figure  15.  Measured  distributions  of  vertical  wind  (left) 
and  vertical  acceleration  of  aircraft  (right): 
shear  produced  (below)  versus  buoyancy  prod¬ 
uced  (above)  turbulence. 


Figure  It.  Measured  spectra  of  vertical  wind  (left)  and  vertical  acceleration  of  air¬ 
craft  (right):  height  dependence. 


Figure  17.  Measured  distributions  of  vertical  wind  (left)  and  vertical  acceleration  of 

aircraft  (right):  height  dependence. 
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Figure  18.  Contributions  of  scale  ranges  to  measured  variances  of  vertical  wind  and 
vertical  aircraft  acceleration.  HP  denotes  total  variance,  B1  contribution  of 
smaff  scales,  B2  contribution  of  intermediate  scales,  TP  of  large  scales. 
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Figure  19.  Comparison  of  3D  LES  simulation  and  aircraft  observations:  vertical  pro¬ 
files  of  vertical  wind  variance. 
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1 4.  Abstract 

-  The  paper  gives  a  brief  overview  of  work  on  atmospheric  turbulence  at  DFVLR.  and  describes 
the  research  tools,  which  include  instrumented  aircraft  and  numerical  models.  Some  results 
of  research  on  turbulence  characteristics  are  given;  these  are  mostly  in  the  convective  boundary 
layer.  The  report  discusses  their  application  to  the  study  of  aircraft  response  and  airframe 
loadings.  ‘ 
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